Majorana Fermions and Non-Abelian Statistics in Three Dimensions 
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We show that three dimensional superconductors, described within a Bogoliubov de Gennes frame- 
work can have zero energy bound states associated with pointlike topological defects. The Majorana 
fermions associated with these modes have non-Abelian exchange statistics, despite the fact that the 
braid group is trivial in three dimensions. This can occur because the defects are associated with 
an orientation that can undergo topologically nontrivial rotations. A feature of three dimensional 
systems is that there are "braidless" operations in which it is possible to manipulate the groundstate 
associated with a set of defects without moving or measuring them. To illustrate these effects we 
analyze specific architectures involving topological insulators and superconductors. 
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A fundamental feature of quantum theory is the quan- 
tum statistics obeyed by identical particles. For ordinary 
particles, Bose and Fermi statistics are the only possibil- 
ities. Emergent excitations in correlated many particle 
systems, however, can exhibit fractional [ij-Ql and non- 
Abelian[5| statistics. The sirnplest non-Abelian excita- 
tions, known as Ising anyons[6j, are Majorana fermion 
states associated with zero energy modes that occur in 
the Bogoliubov de Gennes (BdG) description of a paired 
condensate^. They have been predicted in a variety of 
two dimensional (2D) electronic systems, including the 
V = 5/2 quantum Hall effect [1], chiral p-wave super- 
conductors (SCs) [9*1 and SC-topological insulator (TI) 
structures tlflj. The groundstate of 2A^ Ising anyons 
has a 2^ degeneracy, and when identical particles are 
exchanged the state undergoes a non-Abelian unitary 
transformation llj, [ij] . Recent interest in non-Abelian 



statistics has been heightened by the proposal to use 
these features for topological quantum computation 13 1. 

Fractional and non-Abelian statistics are usually asso- 
ciated with 2D because in 3D performing an exchange 
twice is topologically trivial. In this paper, we show that 
in 3D Majorana fermion states are associated with point- 
like topological defects, and that they obey non-Abelian 
exchange statistics, despite the triviality of braids. Our 
motivation came from the study of 3D SC-TI structures, 
where Majorana fermions arise in a variety of ways, such 
as (i) vortices at SC-TI interfacesllOl, (ii) SC-magnet in- 
terfaces at the edge of a 2D TI[14l4l6| and (iii) band 
inversion domain walls along a SC vortex line. While 
the Majorana fermions in these cases can be identified 
using ID or 2D effective theories, they must occur in a 
more general 3D theory. To unify them we introduce 
a Z2 topological index that locates the zero modes in a 
generic 3D BdG theory. We then study a minimal 8 band 
model in which the defects can be understood as hedge- 
hogs in a three component vector field. Ising non-Abelian 
exchange statistics arise because the hedgehogs have an 
orientation that can undergo nontrivial rotations. We 
will illustrate the intrinsic three dimensionality of the 
Majorana states by considering specific architectures in- 



volving SCs and TIs. A feature in 3D is the existence of 
"braidless" operations, in which the quantum informa- 
tion encoded in the Majorana states can be manipulated 
without moving or measuring 17| them. 

To determine whether a Majorana mode is enclosed in 
a volume V we topologically classify BdG Hamiltonians 
on dV , the 2D surface V . We assume the Hamiltonian 
varies slowly, so we can consider adiabatic changes as 
a function of two parameters r characterizing dV . The 
problem is then to classify particle-hole (PH) symmetric 
BdG Hamiltonians 'H(k, r), where k is defined in a 3D 
Brillioun zone (a torus T'^) and r is defined on a 2-sphere 

. PH symmetry is defined by an anti- unitary operator 
S satisfying = 1 and ■H(k,r) = -S'H(-k, r)E:-i. As- 
suming no other symmetries, this corresponds to class D 
of the general scheme 3 3] • Since V may or may not 
enclose a zero mode, we expect a Z2 classification - a fact 
that can be established using methods of K theory. 

A formula for the topological invariant can be derived 
using a method similar to Qi, Hughes and Zhang's (2^ 
formulation of the invariant characterizing a 3D strong 
TI 21]. We introduce a one parameter deformation 
■H(A, k, r) that adiabatically connects ■H(k, r) at A = to 
a trivial Hamiltonian independent of k and r at A = 1, 
while violating PH symmetry. PH symmetry can then 
be restored by including a mirror image H(A,k, r) = 
-SH(-A,-k,r)S-i for -1 < A < 0. For A = ±1, k, r 
can be replaced by a single point, so the 6 parameter 
space (A, k, r) ^ Y,{T^ x 5^) (S denotes the suspension) 
has no boundary. % defined on this space is characterized 
by its integer valued third Chern character p^. 
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Here F = dA + A /\ A follows from the non-Abelian 
Berry's connection Aij — {ui\d\uj) associated with the 
negative energy eigenstates of T-L. Due to PH symmetry 
the contributions to (P) for A > and A < are equal. 
Moreover, since it is a total derivative, Tr[J^^] = dQ^ 
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(omitting the A's), where the Chern Simons 5-form is(22j 

Qs ^Tt[AA {dAf + {3/2)A^ AdA+ (3/5)^^^] . (2) 

The integral over A > can be then be pushed to the 
boundary A = 0, so that 

A different deformation H can change Ch3[J^], but PH 
symmetry requires the change is an even integer. Like- 
wise, the right hand side of ([3]) is only gauge invariant 
up to an even integer. The parity of ([3]) defines a Z2 
topological invariant. We write it as 

/i = / a;(r) • dA mod 2 (4) 

where the gauge dependent Chern-Simons flux u; = 
{u!^ , uj'^ , uj^) is defined by integrating out k, 

le,,ku^^ir)dx^ Adx'^ = ^(j-^ j^^Q^. (5) 

It is natural to associate /i with the presence of a zero 
mode - a fact that will be checked explicitly below. 

We now introduce a minimal model that leads to an 
appealing physical interpretation for w(r). Since /i is 
based on Ch3[J^], we expect a minimum of 8 bands is 
required. Consider a model parameterized by a three 
component vector field n of the form, 

U = -ijada + TaUair). (6) 

Here ja and Ta {a = 1,2,3) are 8x8 Dirac matrices 
satisfying {rajFf,} = {70, 7b} = 26ab and {Ta,jb} = 0. 
% respects PH symmetry provided SFaS^^ = —Ta and 
S7aS~^ — -fa- For n(r) — no, H has eigenvalues 
E{k) = ±(|k|2 + |nop)i/2^ so that for no 7^ there is 
a gap 2|no|. The seventh Dirac matrix 75 = iJla'Ta-'^a 
is not an allowed mass term because 575S~^ — 75. A 
more general Hamiltonian could also involve products of 
the Dirac matrices, but such Hamiltonians can be homo- 
topically deformed to the form of ([5]) without closing the 
gapJJij. To regularize © at |k| — cx) we include an ad- 
ditional term ejkpFa, so k can be defined on a compact 
Brillouin zone S^. The analysis is simplest for e — !• 0, 
where the low energy properties are isotropic in n. 

H can be physically motivated by considering a BdG 
Hamiltonian describing ordinary and topological insula- 
tors coexisting with superconductivity. The Dirac matri- 
ces are specified by three sets of Pauli matrices: r for PH 
space, (7 for spin and fl for an orbital degree of freedom. 
We identify 7 = t^^^ct, Fi = r^;, F2 = Ty and F3 = t^^j.^, 
along with S = ayTyK. n is then (Ai,A2,to), where 
A = Ai -|- iA2 is a SC order parameter and m is a mass 
describing a band inversion. For A = (jH]) is a doubled 



version of the model for a 3D TI discussed in Ref. 0. For 
e > 0, m > describes a trivial insulator, while m < 
describes a TI with a band inversion near k = 0. An in- 
terface where m changes sign corresponds to the surface 
of a TI, which has gapless surface states. Introducing 
A 7^ to the interface then describes the proximity in- 
duced SC state 0. 

To locate the zero modes, we take n(r) to vary slowly 
with r and evaluate uj{r) using ([TJIS]). This can be done 
by noting that Ji defines a 6 component unit vector given 
by the direction d^(k, r) of {ki, k2, k^, ni + e|kp, 712, 71,3) 
on S^. The deformed Hamiltonian H can be defined by 
adding A75 to H, so that H defines a vector on given 
by = {y/r^d^,X). Ch3[T] is the volume on 5^ 
swept out by d^(A, k, r). Then J Q5 is the volume in the 
"northern hemisphere" of swept out by (k, r) , which 
is confined to the "equator", A = 0. This is then related 
to the area on 5*^ swept out by (k, r) . Performing the 
integral on k for e — >■ gives, 

= ^e'^'n-a,nx a,n, (7) 
on 

where n = n/|n|. Thus, the topological charge that sig- 
nals a zero mode inside V is the parity of the S'^ winding 
number of n on dV. Zero modes are associated with 
hedgehogs in n(r). A simple example of a hedgehog is a 
SC vortex at the interface between a TI and an insula- 
tor. Though the hedgehog topological charge can be any 
integer, an even integer in Q can be unwound by a k 
and r dependent gauge transformation. 

The presence of a zero mode associated with a 
hedgehog can be demonstrated with a simple linear 
model 71a (r) = Mabrb, which has a hedgehog with 
charge sgn(det[M]) at r = 0. This is solved by 
expressing M in terms of its principle axes: M = 
Ofdiag(Mi,M2, M3)02, where Oi and O2 are orthog- 
onal matrices that diagonalize MM'^ and M'^M respec- 
tively. Defining = OiabU, n'a = 02abnb, 7a = ^lab^b 
and Fjj = ^2ab^b, it is straightforward to express "H^ as 
three independent harmonic oscillators, 

=^Afa(277,-f (8) 

a 

where Tia are oscillator quantum numbers and ^a — ^Ta-^a 
are commuting operators. There is a single zero energy 
state with ti^ = and = 1- This zero mode is the 
non degenerate eigenstate with eigenvalue 3 of — 
iY.a,blaOabTb, where O = OiO^ . 

A key feature of the zero mode is its dependence on 
the relative orientation O of the principle axes in r and 
n space. This can lead to a non trivial holonomy when 
M (and hence O) varies. We construct the zero mode by 
starting with |^o) which satisfies X)a *7a-'^al^o) — 3|^'o) 
and then doing a unitary transformation that takes F^ to 
Oab^b- If we parameterize the rotation O with a vector 
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FIG. 1: Inverse image paths depicting the interchange of two 
hedgehogs, as described in the text, (a-e) show a sequence of 
continuous deformations of n(r) that preserve the orientation 
of the hedgehogs. In (f ) a 2n rotation is required to return to 
the original configuration. 



rj specifying the axis and angle |r2| < tt, then 

=e'°''=r"r''"=/4|*o>. (9) 

This gauge satisfies (^'(r2)|d^'(r2)) = for |r2| < tt but is 
not globally defined because = — il)) when 

— TT. This refiects the nontrivial topology of 50(3), 
characterized by the homotopy 7ri(5'0(3)) = Z2. When 
O varies along a nontrivial loop in 5*0(3), the wavefunc- 
tion of the zero mode changes sign. The associated Ma- 
jorana operator 7^ (not to be confused with the Dirac 
matrix 7^) also changes sign. This is one of our central 
results, and it is this fact that allows Ising non-Abelian 
statistics in 3D. A simple example of a nontrivial loop is 
the 2tt rotation that occurs when the SC phase advances 
by 27r. 

Though we derived it with the linear 8 band model, our 
conclusion that there is nontrivial holonomy for defects 
is more general. A formulation based on Q will appear 
elsewhere. In the 8 band model, a general defect his- 
tory is characterized by n(r, t), where r is on a surround- 
ing surface 5^ and t varies on a closed path 5^. These 
are classified by the homotopy of maps 5^ x 5 



^ 

which were first analyzed by Pontrjagin 23| , and have ap- 



peared in other physical contexts [2J]. When the hedge- 
hog's topological charge is ±p the classification is Z2p- 
This is related to the integer Hopf invariant for maps 
— ?> 5^ which applies whenp — 0. Like the Hopf invari- 
ant, it can be understood in terms of the linking of curves 
in 5^ X S-^ 25] . The discussion below resembles Wilczek 
and Zee's|2| analysis of the statistics of skyrmions in the 
2-l-lD nonlinear a model with a Hopf term. 

To study the exchange statistics of the Majorana 
modes we consider the adiabatic evolution of the state 
when they are exchanged. We thus consider a 3-l-lD his- 
tory n(r, t) satisfying n(r, T) — n(r, 0) with hedgehogs at 
{ri{t)} with ri(^2){T) = r2(i)(0). To visualize n(r,i) it is 
useful to consider the inverse image paths in r space that 
map to two specific points on 5^. Such paths begin and 
end on hedgehogs, and a crucial role will be played by 
their linking properties. Fig. 1 depicts four hedgehogs. 



where the top two (positive) hedgehogs are interchanged. 
At the first step in (b) the locations of the hedgehogs 
have been interchanged. Since the inverse image paths 
have been "dragged", n(r) is not the same as it's original 
configuration. Panels (c-f) show a sequence of smooth 
deformations that untangle n(r). The key point is that 
the two paths (which map to different points on 5^) can 
never cross each other. However, a path can cross itself 
and "reconnect", as in (c,d,e). The deformations from 
(a-e) preserve the orientation of the hedgehogs, but leave 
behind a twist. To return n(r) to its original configura- 
tion in (f) requires a 27r rotation of one of the hedgehogs. 
This results in an interchange rule for the Majoranas, 
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71 ^ 72 



72 -71 



(10) 



analogous the rules |lllll2| for braiding vortices in 2D and 
can be represented by T12 = exp [7171 72/4]. Note that in 
(d,e) the twist could have been left on the other side, 
which would have led to T21 = Tj'^2- The two choices for 
T12 correspond to physically distinct interchange trajec- 
tories that generalize the right and left handed braiding 
operations in 2D. 

Performing the same interchange twice leads to a non- 
trivial operation, since T? = 7^7^ changes the sign of 
both 7i and jj . This is natural in 2D because it is a non 
contractable braid. In 3D, however, T^j can be smoothly 
deformed into an operation in which all particles are held 
fixed. Thus there is an operation, specified by a history 
n(r, t), that rotates any pair of stationary hedgehogs by 
27r, and implements the operation 7i7j . The existence 
of such "braidless" operations is a feature of Ising non- 
Abelian statistics in 3D. Although these operations form 
an admittedly limited Abelian subgroup, they nonethe- 
less offer a method for manipulating the quantum infor- 
mation encoded in the Majorana fermions without mov- 
ing or measuring them. 

We now illustrate these effects using specific architec- 
tures involving TIs and SCs. It is easiest to engineer Ma- 
jorana modes using structures involving interfaces or vor- 
tex lines, where u; is confined to lines or planes. Nonethe- 
less, such structures can exhibit intrinsically 3D effects. 
Consider first the structure in Fig. 2(a), which involves 
two disconnected spherical TIs surrounded by a SC and 
connected to each other by a Josephson junction. Sup- 
pose that each sphere has a single ± pair of vortices, so 
that there are 4 Majorana states on the spheres. The in- 
ternal state of the Majorana fermions can be represented 
in a basis of eigenstates of n = 17172 and n' = 17374. For 
an isolated system the parity of n -I- n' is fixed, so the 
system is a single qubit with basis vectors \nn' ~ 00, 11). 
The state can be initialized and measured in this ba- 
sis with a probe that couples to both 71 and 72. Sup- 
pose 73 is adiabatically transported around 71 as shown. 
This is similar to a 2D braid, and it implements T^^g, 
which interchanges |00) and |11). Note, however, that 
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FIG. 2: Architectures demonstrating 3D Majorana states us- 
ing SCs and TIs. (a) Braiding 3 around 1 can be deformed 
into a braidless operation, (b) A geometry for implement- 
ing braidless operations with a Josephson junction device [l5|. 
(c,d) Thin film geometries for interchanging and measuring 
Majorana states, (c) shows a layered SC-TI-SC structure, 
while (d) shows a thin film of a SC weakly doped TL 



this "braid" can be smoothly contracted to zero by slid- 
ing the path around the other side of the sphere. In 
the process, however, the path crosses the junction con- 
necting the spheres, resulting in a In phase slip. Thus 
the braiding operation can be smoothly deformed into a 
"braidless" operation, where 71,73 are held fixed, but the 
phase difference between the two SCs advances by 27r. 

A similar but more feasible version of the braidless 
operation occurs for the Josephson junction structure in 
Fig. 2b, which involves Majorana modes at a SC-magnet 
interface at the edge of a 2D TL As argued in ReL |l5, 
when the phase difference across the junction is advanced 
by 27r, the fermion parity associated with 17172 changes, 
resulting in a fractional Josephson effect. 

Fig. 2c shows a Tl coated on the top and bottom with 
SC films. In a magnetic field, Majorana states occur at 
vortices on both the top and bottom. If the SC films are 
thinner than the penetration depth the field is constant, 
so the top and bottom vortices are independent. If the 
Tl is thin, there will be a weak vertical coupling that 
splits nearby Majorana modes according to n = «7i72 
and n' = 17374. The states \nn') will have slightly differ- 
ent charges, which may allow n, v! to be measured with 
a sensitive charge detector. Suppose the state is initially 
1 00). Interchanging 71, 73 on the top keeping the bot- 
tom fixed leads to the entangled state (|00) -|- |ll))/-\/2. 
A variant on this geometry (Fig. 2d) is a thin film of a 
bulk SC weakly doped Tl, in which the surface states ac- 
quire SC similar to the proximity induced state. In this 
case, the interchange of 71, 73 involves a reconnection of 
the vortex lines (which in principle have a finite energy 
gap) connecting them. 

Receiitly, SC has been observed in Cua;Bi2Se3 for 
X ~ .15[26j. It will be interesting to determine whether 
this material is in the weakly doped regime, with Majo- 
rana modes at the ends of vortex lines, or a more con- 



ventional SC, which could be used in Fig. 2c. There 
are certainly technical challenges associated with manip- 
ulating the vortices and measuring their charge state. 
Nonetheless, we hope that the prospect of detecting 3D 
non-Abelian statistics in such a system will provide mo- 
tivation for further exploration. 

We thank Liang Fu for insightful discussions and Bryan 
Chen and Randy Kamien for introducing us to the Pon- 
trjagin invariant. This work was supported by NSF grant 
0906175. 
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